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We study the process of production of ultracold molecules from ultracold atoms using a sinu- 
soidally oscillating magnetic field modulation. When the magnetic field is resonant roughly with 
the molecular binding energy, Shapiro-like resonances are observed. Their resonance profiles are well 
fitted by the Lorentzian functions. The line widths depend on both the amplitude and the duration 
of the applied modulations, and are found to be dramatically broadened by thermal dephasing effect. 
The resonance centers shift due to both many-body effect and finite temperature effect. Our theory 
is consistent with recent experiment (S. T. Thompson, E. Hodby, and C. E. Wieman, Phys. Rev. 
Lett. 95, 190404 (2005)). Our model predicts a 1/3 ceiling for the molecular production yield in 
uncondensed ultracold atomic clouds for a long coupling time, while for the condensed atoms the 
optimal conversion yield could be beyond the limit. 

PACS numbers: 03.75.Nt, 34.50.-s, 36.90.+f 



I. INTRODUCTION 



Shapiro resonance is one of the most remarkable prop- 
erties of the superconducting device, in which, two 
weakly coupled superconductors are subject to a voltage 
difference that is the sum of a dc component V and a pe- 
riodic signal V m s'm(ft). A continuous range of nonzero 
dc currents are possible if V = j-kf, where 2e is the 
Cooper pair charge, h is the reduced Planck constant, 
and k is an integer [jj, Q. The Shapiro resonance provides 
a method to measure the constant of nature 2e/h with 
such precision and universality @,IH that, since 1972, the 
reversed view has been adopted whereby 2e/H is assumed 
to be known and the above Shapiro resonance is used to 
define a standard unit of voltage 0, [H, @. 

Essentially, the Shapiro resonance is a specific phe- 
nomenon emerged when the frequency of the external 
field is commensurate with the intrinsic frequency of sys- 
tem. Recently, it has received renewed interests and in- 
vestigations in the Bose- Einstein condensates(BEC)@,[l, 
9]. For example, in BEC Josephson junction, the dc 
value of the drift current shows up as resonant spikes 0- 
Under experimentally accessible conditions there ex- 
ist well-developed half- integer Shapiro- like resonances Q. 
Shapiro effect also allows precise measurements in atomic 
BECs. The ac-driven atomic Josephson devices can be 
used to define a standard of chemical potential^. 

In the present paper, we extend to investigate the 
Shapiro resonance effects in ultracold molecule produc- 
tion. The conversion of ultracold atoms to ultracold 
molecules by time varying magnetic fields in the vicinity 
of a Feshbach resonance is currently a topic of much ex- 
perimental and theoretical interest. This particular con- 
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version process lends itself well to the formation of molec- 
ular Bose-Efnstein condensates (BECs)[H [H 0J, El 
and atom- molecule superpositions (hi] . These Feshbach 
molecules and their creation process are also impor- 
tant for understanding ultracold fcrmionic systems in the 
BCS-BEC crossover regime because they are closely re- 
lated to the pairing mechanism in a fermionic superfluid 
that occurs near a Feshbach resonance [H, [HI, [l?], Hcj |. 
We study the process of production of ultracold molecules 
from ultracold atoms using a sinusoidally oscillating mag- 
netic field modulation. The advantage of this method 
is that it greatly reduces the heating the cloud expe- 
riences in the conversion process because the conver- 
sion occurs far from the center of the Feshbach reso- 
nance. In recent experiments, this technique has been 
applied and was shown could produce molecules from 
atoms more emcientlypjl and measure the binding en- 
ergy of Feshbach molecule precisely. However, the un- 
derlying mechanism is not fully understood, the com- 
plexity arises from the many-body problem and time- 
dependent field involved. The many experimental ob- 
servations can not be accounted for by the existing the- 
oretical model. For example, the Poisson distribution 
predicted by theory 20] is departure from the observed 
Lorentzian-likc resonance profiles and how the finite- 
temperature effect influences the resonance profiles is still 
an open problem. In the present paper, we exploit a 
microscopic two-channel model to investigate thoroughly 
the mechanism underlying the Shapiro resonance phe- 
nomenon in the atom-molecule conversion. With quan- 
titatively considering the thermal dephasing effect in the 
uncondensed atom clouds our model could account for 
the most experimental observations. Our theory also sug- 
gests some interesting predictions for future's experimen- 
tal test. 

The plan of this paper is as follows. In Sec. II we 
present our model and make a thorough analysis on the 
Shapiro resonance. In Sec. Ill, with the inclusion of the 
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dephasing effect in our model, we apply our theory to 
explain the recent experiment. In Sec. IV, we extend to 
discuss the Shapiro resonance for the case of condensed 
atoms with emphasizing on the influence from the inter- 
action between coherent particles. Sec. V is our conclu- 
sion. 
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II. SHAPIRO-LIKE RESONANCE IN 
ATOM-MOLECULE CONVERSION FROM A 
TWO-CHANNEL PERSPECTIVE 
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A. Model 

Ignoring the two- and three-body atomic decay and 
collisional molecular decay, we exploit the following two- 
channel microscopic model to describe the dynamics of 
converting atoms to molecules in the bosonic system, 

H = (e a - n) a) a + (e b + v{t) - 2/j) 

+-^=(a^a r b + Paaj . (1) 

Here a (a + ) and b(b + ) are bose annihilation (creation) 
operators of atoms and molecules, respectively. The to- 
tal number of particles N = at a + 2w b is a conserved 
constant. The atomic and molecular kinetic energies are 
given by e a and £(,, fi is the chemical potential, g gov- 
erns the atom-molecule coupling strength, V denotes the 
quantization volume of trapped particles and therefore 
n = N/V is the mean density of initial bosonic atoms. 
Where v(t) represents the binding energy of diatomic 
molecules which depends on the external field, expressed 
approximately as[2l| . 



u(t) = - 



(2) 



where ro is the effective range of the van der Waals po- 
tential, to is the mass of a bosonic atom, and a e ff de- 
notes the effective scattering length driven by external 
magnetic field, 
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where ab g is the background scattering length, Bq is the 
Feshbach resonance position, AB is the width of the 
resonance defined through the relation with the atom- 
molecule coupling term AB = mg /Airh \aj, g AiJ,\, where 
Afi is the difference in magnetic moment between the 
closed channel and the open channel state. We focus on 
the situation that the selected external field B ex is mod- 
ulated sinusoidally with small amplitude B mo( i and large 
frequency u> near a Feshbach resonance, i.e., 



B(t) = B ex + B mod sm(u)t). 



(4) 



Since B mo d -C B ex , the binding energy can be expanded 
into series to the first order of B mnr i 



v(t) = v e +v m sin(wi), 



(5) 



B. Shapiro-like resonance 

We introduce the angular momentum operators to in- 
vestigate the dynamics of this system[22[, 
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Then the Hamiltonian can be written as 
N 



H 



[iyo + v m sia(u>t))L z + V2i]L x ^ , (14) 



where vq = v e + £b — 2e a is the energy difference between 
atoms and molecules, and parameter rj = 2gy/n denotes 
the coupling strength. Then the Heisenberg equations of 
motion are 
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Since N is large for the current experiments and all the 
commutators vanish in the limit of N — > oo, it is ap- 
propriate to take L x , L yi and L z as three real numbers 
u, v, w, respectively. Then we get the mean-field Heisen- 
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From the Heisenberg equations (jl 511 17[) , we get 
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FIG. 1: Time-averaged population imbalance —{L z )t for the 
driven system with different numbers of particles, tilt vo/r/ = 
5, and scaled driving amplitude v m = 1. 
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In order to get the time-averaged value of the conversion 
varied with different external field, we characterize each 
quantum trajectory by its time-averaged imbalance 
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employing the averaging interval At 3> fi/vo- Initially, 
all particles are atoms. FigQ] shows the results of such 
calculations by numerical solving the Heisenberg equa- 
tions (JTSHT]) for N — 2, 20 under periodic modulation 
with fixed scaled amplitude v rn jva = 0.2 and frequen- 
cies u> ranging from to 1.25^o- The solution of the 
mean-field equations 



!)18H20p is also presented. There are 
several clear spikes which indicate the Shapiro-like res- 
onance in atom-molecule conversion driven by external 
magnetic field. 

These spikes indicate that the frequency of the modu- 
lated field is commensurate with the intrinsic frequencies 
of the atom-molecule conversion system in the absence 
of the periodic modulation. Now we analyze the intrin- 
sic frequency. For N = 2, using Fock state as basis, the 
commutators (|lltil3[) becomes 
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Then the intrinsic fre quency i s readily obtained from the 
above equation as ^/v 1 + rj 2 /h. So the center of reso- 
nance is expected to be ^Jv^+rf / (ftio) — p/q with p, q 
are integers. In our case, the resonances corresponding to 
p/q = 1,2,3 are more prominent. With N increasing, we 
find that the resonance center shifts to right due to the 
many-body effect. We can obtain the intrinsic frequency 
in the mean-field limit, i.e., N — > oo. From the mean-field 
equations (|18H20p we readily obtain, 



0. 



(25) 



Initially all particles are in atom states, i.e., w = 
— 1. Approximately substituting it into the above 
equation we obtain the explicit expression of the fre- 
quency y v 2 + Arj 2 /H. It implies that due to the 
many-body effect [231] . the resonance centers shifts to 
\J + A-q 2 / (ftw) = 1, 2, 3, .... The above theoretical anal- 
ysis agree with our numerical results. 



C. Phase space at the Shapiro-like resonance 

The Shapiro resonance phenomenon can be demon- 
strated intuitively by the trajectories in phase space 
of the system. Notice that the constraint u 2 + v 2 — 
-(«; — l) 2 (w + 1) and introducing the canonical vari- 
able s = w, 9 — arctan(i>/u) denoting the population 
imbalance and the relative phase between atoms and 
molecules, the mean-field Heisenberg equations can be 
replaced by a classical Hamiltonian of the form 
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the canonical equations of motions are 
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FIG. 2: Poincare section of the classical Hamiltonian (|26l 
with tilt vo/r) = 24, scaled driving amplitude v m = Q.2vq 
and modulation frequency (a) uj/vq = 0.95 (off-resonance) 
(b) lu/vq — 1 (near resonance). 



wc obtain the new Hamiltonian 
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The secular evolution of S and can be evaluated from 
the time-averaged Hamiltonian, 



(JC(s,e)) T = i £ ic(s 7 Q)dt, 



(33) 



with T = 2tt/lu. 

Now, we consider the 1:1 resonance case that hw sa uq, 



(JC(S,e)) T « Ji( — )?V(1 - S) 2 (l + 5) sine. (34) 

where Ji(x) is the first kind Bessel function. The phase 
graph of the periodic averaged Hamiltonian system re- 
flects the Poincare section of Hamiltonian system (|26[). 
as shown in Fig [2] For the case of off-resonance, the inte- 
gral in Eq. (|33[) approximates to zero, it implies that the 
time-averaged s varies a little in time while the variable 
6 increases with time almost linearly. The corresponding 
Poincare section is shown in FigJSJb). ft is shown that 
the phase space at the transition change dramatically. 



D. Intrinsic resonance width: Arnold tongues 

For the initial condition s = w — 1, whether its tra- 
jectory falls into resonance regime can be judged from 
following resonance condition AH > 2vq/H, where A7Y is 



FIG. 3: Arnold tongues for resonance modes (1 : 1),(2 
1), (3 : 1), (4 : 1), (5 : l)(from right to left) with v /v = 5. 



the difference between the maximum and the minimum 
value of H in the time interval At ^> H/vq- For different 
v m and uj, we obtain the regions in the two-dimensional 
parameter space where the resonance emerges. These re- 
gions are named as Arnold tongues (2~il|. In order to draw 
out the Arnold tongues in parameter space, the main nu- 
merical tool used in this work is the winding number W 



W 



lim 

t — >oo 



6(t) - 9(0) 



t 



(35) 



with initial conditions (0(0), s(0)). If the ratio W/w is ra- 
tional, i.e., W/u> = q/p, here q andp are natural number, 
(6(t),s(t)) is a resonant solution of (q : p) type, i.e., 

(6(t + pT),s(t + P T)) = (6(t), s(t)) + (2th?, 0), (36) 

which means the system runs q times in time pT interval. 
In our system, only (q : 1) type is significant. In Fig[3] 
we show the first five resonance regions with Vo/w = 5. 
The width of the resonance regions is broadened as the 
modulation amplitude v m increases. 



III. COMPARISON WITH EXPERIMENT 

In the above discussion we use the single-mode model 
to discuss the conversion between the condensed atoms 
and molecules. This is an approximation because in prac- 
tical experiments the atoms are not condensed and other 
modes will be coupled. We use this approximation under 
the condition that the energy distribution of the ther- 
mal particles (characterized by ksT, ks is the Boltz- 
man constant and T is the temperature) is much smaller 
than the effective Feshbach resonance width ff^/njl^]. In 
such cases, each 'energy band' of the thermal particles 
can be approximately denoted by one energy level, as 
schematically plotted by Fig0J Initially, the particles 
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FIG. 4: Schematic of the swept magnetic field in experiment. 
See text for details. 



on one level have a definite phase and the phase differ- 
ence between two levels is well defined. However, as the 
magnetic field sweeps and then the sinusoidal magnetic 
field applied, particles will acquire additional phases that 
are proportional to their individual energy and evolution 
time. The varied particles in one level could acquire dif- 
ferent phases because they have different energies. This 
define a 'dephasing rate' 7 = fcnT 7 '(27rfe)[2a|. To to com- 
pare with experiment, we need to include the dephasing 
effect into our model. Modeling dephasing by fully in- 
clude the quantum effects requires sophisticated theoret- 
ical studies. The standard approaches of quantum optics 
for open systems involve quantum kinetic master equa- 
tions. Here, we adopt the simple mean-field treatment 
in our model. From the mean-field viewpoint, the deco- 
herence term introduces a 7 transversal relaxation term 
into the mean- field equations of motion [25l|. 



dt 
d 

—1 

dt 



dt 



-jiy{t) + £h - 2e a )v ~ ju, 



= -(v(t) + e b - 2e a )u 

+ f^ ( »-i)(» + i 



(37) 

(38) 
(39) 



The imbalance of atom-pairs and molecules w is varied 
in the range of [— 1, 1] with the lower limit corresponding 
to a pure atomic gas and w = 1 for a pure molecular gas. 
What we concern is that after the conversion process how 
many atomic pairs are converted to molecule. We use Wf 
to denote the value of w when the magnetic field sweep 
back. The molecular conversion efficiency can be read 
from the variable Wf as T = (1 + Wf)/2. 

Now we apply our theory to the experiment of 85 Rb 
by Ref.[19]. The atoms are held in a purely magnetic 
trap at a bias field of B r . After evaporative cooling, the 



FIG. 5: Conversion efficiency of atoms converted to molecules 
as a function of modulation frequency for three different cou- 
pling times. For a fixed coupling time, the curve can be fitted 
by a Lorentzian distribution T = To + ^-J^yi+A' 2 ' e '6-> 
for t c = 38ms, the fitting parameters are Fo = 0.06, uj c = 6.2, 
A = 0.6, A = 0.23. In the subfigure, by fitting the linewidth 
versus conversion data to a straight line we find the zero con- 
version limit to be 0.3kHz. 



magnetic field B is linearly swept to a selected value at 
B ex and then apply a sinusoidal magnetic field pulse with 
peak-to-peak amplitude B mo d and modulation frequency 
lo for a duration of coupling time. The swept magnetic 
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FIG. 6: Conversion efficiency of atoms converted to molecules 
under a periodic modulation with amplitude B mo d = 0.13G 
and frequency u) = 6.25kHz with respect to coupling time 
for different temperatures. The conversion of ultracold atoms 
to molecules increases with coupling time until it becomes 
saturated at 1/3. 
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field can be expressed as 



B r -at < t < t 

B = { B ex + B mod sin(wt) t < t < t + t c 

B ex + at t + t c < t < 2t Q + t c . 

(40) 

Where B r = 162G, B ex = 156. 5G, B rnod = 0.13G, w 
is ranging from 2kHz to 9kHz, to is the linear sweep 
time, a = (B r — B ex )/to is the linear sweep rate, and 
t c is the coupling time. The sketch curve is shown in 
FigJU For the thermal cloud, with temperature T, one 
molecule have 5 degrees of freedom while two atoms have 
6 degrees of freedom, according to the equipartition theo- 
rem, we have (2e a — £&) « fcsT/2. The scaled parameters 
in Eq. §7Mfr are 



!/(<) 



h 2 



(B - B ) 2 



and 



rj = 2j4irh 2 \a bg Au\ABn/7 



(41) 



(42) 



The experimental parameters are a bg — — 443ao, ro — 
185a [|j|, AB = 10.71G, B = 155G, A/x = 1.2 x 
10~Vb, the temperature T = 20nK, density n = 
10 n cm -3 , here ao and ug are Bohr radius and Bohr 
magneton, respectively. The difference of magnetic mo- 
ment A/j, is extracted from the experimental data 27 1. 
Under this condition, the ratio between the energy differ- 
ence and the energy distribution of thermal particles, i.e., 
Ve/ksT, is estimated to be around 15 that is much larger 
than one. The above analysis validates our single mode 
approximation. Fig|5]shows the conversion efficiency as a 
function of modulation frequency for three different cou- 
pling times. The resonance line width are broadened by 
the dephasing term. There is a clear Lorentzian distribu- 
tion resonance at frequency about 6.25kHz, close to ex- 
periment. Except the fundamental frequency resonance 
at oj = 6.25kHz, there is also a weakly (2:1) mode 
resonance at about u = 3.1kHz, while it has not been 
observed in experiment. Our linewidth is approximately 
0.3kHz at zero conversion limit, as shown in the subfig- 
ure. In the experiment, it is about 0.2kHz. 

In our calculation, we find that, in the three stages of 
magnetic field change expressed by Eq. (|4"0")) . the linear 
process contributes little to the atom-molecule conver- 
sion. This is because oscillation center B ex is still far 
away from the Feshbach resonance center. The atom- 
molecule conversion mainly occurs in the process of ap- 
plying the sinusoidal magnetic field, where v(t) can be ex- 
pressed as vo + Vm sin(wt). Therefore, the above observed 
resonance phenomenon corresponds to the Shapiro reso- 
nance discussed in the last section, while the linewidth is 
dramatically broadened by the thermal dephasing effect. 

In Fig [HI we show the conversion efficiency with re- 
spect to coupling time. For temperature T — 20nK, den- 
sity n = 10 n cm -3 , our results are close to experimental 




2 4 6 

Coupling Time (ms) 

FIG. 7: The conversion efficiency from condensed atoms to 
molecules under a periodic modulation with fixed amplitude 
B mo d = 0.5G and different frequencies . The density for the 
condensed atoms is n — 10 12 cm -3 . The dark triangle marks 
the experimental observation in Ref . [l9| . 



data. We also show the cases of different temperatures 
by considering the isobaric condition, i.e., nT — const.. 
The above calculation shows that increasing the temper- 
ature will lessen the molecular production because the 
dephasing term is proportional to temperature. On the 
other aspesct, the conversion efficiency decreases with 
the increasing of temperature. For different tempera- 
ture, a common feature is that the conversion efficiency 
increases with coupling time until it becomes saturated 
at 1/3. This can be explained from investigating Eq.(37- 
39), where u = v = 0,w = —1/3 is the fixed point in the 
absence of the dephasing term. In the presence of the 
dephasing effect, the system is expected to relax into the 
fixed point in the case of long-term coupling. 



IV. DYNAMICS OF ATOM-MOLECULE 
CONVERSION FOR THE CONDENSED ATOMS 



In this section, we extend our discussion to the case 
of condensed atoms. For pure BEC atoms, the single 
mode approximation is valid and dephasing effect can be 
ignored while the interaction between the coherent atoms 
become significant. After ignoring the kinetic energies of 
particles, the Hamiltonian can be written as 



H = v(0b - — a+S+aa + -j= (tftfb + Pa&j .(43) 



where U — ATrh 2 \ab g \/m denotes the nonlinear interac- 
tion. In mean-field limit, we can derive the Heisenberg 
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CO (kHz) 

FIG. 8: The conversion efficiency from condensed atoms to 
molecules with respect to modulation frequency for different 
nonlinear interaction U. The coupling time is 1.6ms. The 
other parameters are the same as in Fig0 



0.14. One sees that there is a Rabi oscillation which can 
reach a high conversion efficiency. In experiment, they 
observe 55% conversion for coupling time 1.6ms jl9|. The 
observation marked by a dark triangle in FigJT] is consis- 
tent with our result. Figj8] presents the conversion effi- 
ciency with respect to modulation frequency for different 
nonlinear interaction U . The coupling time for this cal- 
culation is 1.6ms. The other parameters are the same as 
in FigUl 

In Fig. 8, we observe some oscillations except for the 
main resonance peaks and find that the maximum con- 
version efficiency could be far beyond the limit 1/3. Be- 
cause the resonance center can be still approximated by 
OJ = \Jv~l + 4rj 2 /h. It implicitly depends on the density 
through parameter 77. In the experiment, the density of 
condensed atoms is about ten times larger than that of 
thermal cloud, therefore the resonance center shift to the 
right-hand compared to Fig. 5 of the thermal atomic cloud 
case. Moreover, we find that the interaction between the 
coherent atoms will lead to the shift of the resonance 
profile as clearly shown in Fig. 8. 



equations of motion 



dt 
d 

—1 

dt 



dt 



u(t) 2Un 
--— 1; z~ v V- ~ w )' 



2Un 
H — z-u(l - w), 
n 

n 



l)(w 



(44) 

(45) 
(46) 



FigO presents the conversion efficiency under a pe- 
riodic modulation with fixed amplitude B moc i = 0.5G 
and different frequencies by numerical solving Eq. (|44H46| . 
The density for the condensed atoms is n = 10 12 cm -3 
and thereby the scaled nonlinear interaction is Un/rj — 



V. CONCLUSIONS 

In conclusion, we have investigated thoroughly the 
mechanism underlying the Shapiro resonance phe- 
nomenon in the atom-molecule conversion with exploit- 
ing a microscopic two-channel model. With the inclusion 
of the thermal dephasing effect in the uncondensed atom 
clouds our model could account for the most experimen- 
tal observations. We also extend our discussions to the 
case of condensed atoms. Our theory have some interest- 
ing predictions waiting for future's experimental test. 
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